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1. Introduction 

Thin-film ferromagnetic materials have played a central role in information storage 
technologies for many years [1-3]. In this context, much attention has been devoted 
to studying these materials by both the physics and mathematics communities 
[1,4]. Magnetic storage media make use of magnetic domains —regions of uniform 
magnetization separated by thin transition layers called domain walls —to represent 
bits of information. Typically, thin films possess one distinguished in-plane direction 
along which the magnetization prefers to align; this direction is referred to as an easy 
axis, and materials possessing a single easy axis as uniaxial. In such a material, this 
results in two distinct optimal domain orientations, which may be used to store binary 
information. More recently, more preference has been given to perpendicnlar materials 
— magnetic storage materials in which the easy axis is in the ont-of-plane direction [2,5]. 
Nevertheless, thin film materials with in-plane anisotropy continue to be important 
for many applications, such as magnetoresistive random access memory [6-10] and 
spintronics [11]. 
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In this article, we study one-dimensional domain walls in thin ferromagnetic dims in 
which the magnetization is strongly penalized from pointing out of the film plane, with 
two orthogonal in-plane easy axes (and thus four optimal magnetization directions). 
We refer to these materials as ultrathin fourfold films. Such behaviour is, for example, 
experimentally realized in very thin (3-19 monolayers thick) films of epitaxial cobalt and 
gives rise to unusual magnetic domain morphologies in which the magnetization vector 
has a tendency to rotate by integer multiplies of 90° in the him plane [12,13]. 

In order to understand the magnetization behaviour of thin-hlm materials from a 
theoretical point of view, one would like to study the formation of domain patterns, 
and the structure of the domain walls which connect them. We start from the Landau- 
Lifshitz-Gilbert (LLG) equation for the dynamics of the magnetization vector M(x, t), 
which is of hxed length |M(x, t)| = Mg and dehned on a spatial domain G C 
representing the ferromagnetic body under study. The LLG equation for M reads 

dM 7 

^ = --^(MxH + «MxMxH), h-VM|an = 0, (1) 

at 1 

where 7 is the gyromagnetic ratio, a is the dimensionless Gilbert damping parameter, 
and H is the effective magnetic held. This held is obtained via H = —6E/6'h/L, where 
T'(M) is the micromagnetic energy functional: 




+ 



dh 

V-M(r)V-M(rO ^ 3 ^ ^ 3 ^, 


( 2 ) 


Here, A is the exchange constant; K is a crystalline anisotropy constant, with <I)(M) a 
scalar function describing the anisotropy; and G is an external applied magnetic 
held. The terms in the energy may be understood as follows. The hrst term is the 
exchange energy, which penalizes spatial variations of M; the second is the anisotropy 
energy, which describes the preferred directions for M within a material; the third is 
the Zeeman energy, which prefers M to align with the external held; and the fourth (in 
which M is extended by zero outside of G and the derivatives are understood in the 
distributional sense) is the nonlocal magnetostatic energy, which prefers to minimize 
the distributional divergence of M. See e.g. [1] for further details of the micromagnetic 
model. 

We can consider observable static domain patterns as stationary (i.e. with ^ = 0) 
solutions of (1), and note that, in this case, equation (1) coincides with the Euler- 
Lagrange equation for E incorporating the pointwise constraint |M| = Mg. This will 
enable us to employ variational techniques to characterize solutions. 

We consider a reduction of the micromagnetic theory appropriate for very thin 
hlms. Many such reductions have been presented before [14-18], corresponding to a 
variety of regimes of the physical parameters in the energy (2). In order to understand 
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the parameter regime we study here, we introduce the following quantities: 


f—) 

V47rM2 J 


1/2 




( 3 ) 


respectively called the exchange length, Bloch wall width, and quality factor. In 
extended dims, we may take the spatial domain hi = x (0,d), where d is the 
him thickness. The physical regime we consider is then characterized by the scalings 
d i L, with LdjO^ ~ 1. This regime (ultrathin, moderately soft him) is relevant for 
a variety of materials [19]. In this regime, one may introduce the dimensionless thin-film 
parameter 

AnM^d Ld d 

which characterizes the strength of the magnetostatic interaction relative to both 
exchange and anisotropy. 

One may then formally derive a reduced LLG equation from (1) by considering the 
limit Q —)■ 0 and d —)■ 0 together with n = 0(1) hxed [20,21]. Letting m = M/M^ and 
assuming m = (mi, m 2 , 0), i.e., that m lies entirely within the him plane, after suitable 
rescalings one hnds the ehective overdamped equation 


cim 

—— = — m X m X h. 

dt 


( 5 ) 


Here, h is the ehective held obtained now as h = —58/dm, where 8 is the reduced 
thin-hlm energy 


£^(m) = ^ / (|Vm|^ + <I)(m)) ^ j [ 

4 Jr 2 «7r J^2 Jr; 


V • m(r)V • m(r'; 2 12 / 


r — r 


drdr', (6) 


and now m : ^ is the in-plane magnetization direction held. In what follows, we 

consider the case of 

<I)(m) = (m ■ ei)^(m ■ es)^, (7) 

corresponding to fourfold anisotropy. This type of magnetocrystalline anisotropy is 
very common for ultrathin epitaxial hlms [19]. Ultrathin hlms with this type of 
anisotropy has been proposed for applications to multi-level magnetoresistive random 
access memories [22,23] and could be of interest to domain wall based devices [10,24]. 

In uniaxial thin hlms, the behaviour of 180° domain walls has been extensively 
studied. For simple one-dimensional prohles connecting the two optimal directions, there 
are two possibilities: the Bloch wall, in which the magnetization transitions between the 
optimal domains by rotating out of the him plane, and the Neel wall, in which the 
rotation occurs entirely within the plane. Which wall type is energetically preferred 
depends essentially on how severe the penalty for rotating out of plane is, which in turn 
depends on the him thickness. In the ultrathin regime we consider in this article, this 
penalty is strong enough to simply forbid any out of plane component of m, so that the 
Neel wall prohle is the only choice. 
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Neel walls have been studied for many years, with a degree of controversy (see 
e.g. [1,25]). More recent micromagnetic studies have led to a good present understanding 
of the Neel wall’s internal structure [1,4,26-31], the main features of which (sharp inner 
core with slowly decaying tails) have been verihed experimentally [32-34], 

Rigorous mathematical studies of the Neel walls began with the work of Garcia- 
Cervera, who studied, both analytically and numerically, the one-dimensional 
variational problem obtained from the full micromagnetic energy by restricting to 
prohles which depend only on one spatial variable [26,27]. The same functional was 
studied by Melcher, who restricted the admissible magnetization conhgurations to 
those constrained to the him plane, and established symmetry and monotonicity of 
minimizers connecting the two optimal directions [28]. Uniqueness of the Neel wall 
prohle and its linearized stability with respect to one-dimensional perturbations was 
treated by Capella, Otto and Melcher [29]. Stability of geometrically constrained 
Neel walls with respect to large two-dimensional perturbations has been demonstrated 
asymptotically in [35]. Most recently a comprehensive study of Neel walls under the 
inhuence of applied magnetic helds was undertaken by Chermisi and Muratov [30]. 
They proved existence, uniqueness, strict monotonicity and smoothness of the wall 
prohle along with estimates for its asymptotic decay. 

To summarize, in ultrathin uniaxial hlms the magnetization is ehectively 
constrained to lie completely in the him plane, and one encounters 180° Neel walls 
as the optimal transition layer prohles connecting the two uniform states. These are 
now well understood. Beyond that, it is possible to observe stable winding domain 
walls, in which the magnetization makes a number of full 360° rotations (most often 
just one, though more are possible) [36-38]. This type of domain walls has received 
recent theoretical attention in [39,40]. A reproducible way to inject 360° walls into 
ferromagnetic nanowires and successful manipulation of such domain walls were recently 
demonstrated experimentally in [41,42]. 

In fourfold hlms with in-plane magnetizations, we can make the following analogies 
with the uniaxial case. In fourfold materials, the 90°-walls are expected to exist as 
optimal prohles connecting two adjacent minima of the potential <I> (e.g. -|-ei and - 1 - 62 ). 
This is analogous to the 180° Neel walls in uniaxial materials. For a 180°-wall in a 
fourfold material, the magnetization has to connect two nonadjacent minima of <F while 
passing directly through a third somewhere in between (i.e. connect -|-ei and —ei, while 
passing through - 1 - 62 ). Moreover, this should occur without the wall simply splitting 
into two separate 90°-walls. This is analogous to the 360°-walls in uniaxial materials. 

In this article we extend the methods contained in previous work concerning 180° 
and 360° domain walls in uniaxial materials to the setting of fourfold materials, and 
prove existence results for both 90° and 180° walls in these materials. These walls, 
despite some apparent analogies with those found in uniaxial hlms, have not been 
previously investigated theoretically. 
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1.1. Reduced model for one-dimensional domain walls 


Since stationary solutions of (5) coincide with critical points of ( 6 ), in order to study 
stationary one-dimensional domain wall profiles, we now seek to derive a ID variational 
problem from ( 6 ) which is appropriate to capture such prohles via minimization. 

In what follows we explicitly restrict to stationary prohles, m(xi, a; 2 , t) = m(a;i, a; 2 ). 
It is convenient to introduce the in-plane magnetization angle 6 ^ : —)■ M via 


m = —ei sin 6 -\- e 2 cos 6. 


( 8 ) 


We now assume a one-dimensional prohle 6 *(a;i,X 2 ) = 0{f) varying only along the 
direction = ei cos(3 -\- 62 sin/9; we refer to the angle (3 as the wall orientation. With 
these assumptions, the LLG equation (5) for a stationary ID prohle 6{x) reduces to 

Q =-9^^ +-smW+ ^cos{0 - (3) sin(6'-/3), (9) 


/ 2 \ V2 

where (—^) is the negative ID half-Laplacian (a linear operator from H^{ 
modulo additive constants, to L^(M) whose Fourier symbol is |/c|). Equation (9) is 
also the Euler-Lagrange equation corresponding to the energy 

( 10 ) 


^/3W = 2 da;+ -||sin(0-/3) 11 ^ 1/2 

where we introduced the homogeneous (semi-)norm [43]: 




1 



{u{x) - u{y)f 
{x - yY 


da; dy. 


( 11 ) 


It is not too difficult to see that this energy corresponds to the energy ( 6 ) of a ID prohle 
per unit width in the transverse direction. 

This model forms the basis of the rest of this article. It is necessary to specialize 
it further to individually examine the two types of wall we study: 90° and 180° walls. 
To state our main results, we need to introduce the following general admissible class 
of functions for a given a G M: 


Ao^ = {ee HlyW) -.e-y^e H\R)}, 


where r]a £ C'°°(M) is a hxed function which satishes 


r]a{x) 


a for a; G (—00, — 1 ), 
0 for X G (1, -|-oo). 


It is easy to see that the dehnition of the admissible class Aa does not depend on the 
specihc choice of [30]. Also, by Morrey’s theorem we may always assume that if 
6 G Aa, then 6 G C(M) fl L°°(R) and that lima-^+oo 6*(a;) = 0 and lima,^_oo 6'(x) = a. 

In the following, we will look for minimizers of the one-dimensional domain wall 
energy in the admissible classes Aa with a = 7r/2 and a = vr to study 90° and 180° 
walls, respectively. 
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2. Main Results 

We are now in a position to state the main resnlts of this work. The hrst result below 
concerns 90°-walls, and provides existence of these as energy minimizing conhgurations. 
These prohles only exist for a particular orientation (3 = —tt/A (modulo 7 r /2 rotations). 
Furthermore, we are able to extract further information on the prohles including 
uniqueness, smoothness, and strict monotonicity. 

Theorem 1 (90°-walls: existence, uniqueness, regularity and strict monotonicity). For 
(3 = —71/A and each u > 0, there exists a minimizer of the energy over the 

admissible class At,/ 2 - The minimizer is unique (up to translations), strictly decreasing 
with range equal to ( 0 , 7 r/ 2 ), and is a smooth solution of (9) satisfying the limit 
conditions 

lim 9{x) = 0, lim 6{x) = 7 r/ 2 . (12) 

x^-\-oo X—)•—oo 

Moreover, if 9^^'^ : M —)■ ( 0 , 7 r/ 2 ) is the minimizer of S_.„/4{9) over At^ 12 satisfying 
6'^°)(0) = 7 r/ 4 , then 9^^\x) = 7 r /2 — 9^^\—x). 

The choice of the admissible class At^i 2 serves to enforce the asymptotic behavior in (12). 
We note that for other wall orientations (3 the wall would carry a net line “charge” and, 
hence, the last term in the one-dimensional wall energy ( 10 ) would always be inhnite 
on An/2- 

The next result is similar to the hrst, but concerning 180°-walls. Again the 
orientation of the walls is restricted, this time to /3 = 0 (modulo 7 r /2 rotations). Many of 
the properties of 90° walls follow here, though uniqueness of the prohle is not presently 
clear. 

Theorem 2 (180°-walls: existence, regularity and strict monotonicity). For (3 = 0 and 
each iz > 0, there exists a minimizer of the energy £p{9) over the admissible class An- 
The minimizer is strictly decreasing with range equal to (0,7r), and is a smooth solution 
of (9) satisfying the limit conditions 

lim 9{x) = 0, lim 9{x) = tt. (13) 

X^ + OD X^ —CXD 

Moreover, if9^^'> : M —)■ (0, vr) is the minimizer of£o{9) over An satisfying 6*’^°)(0) = 7t/2, 
then 9^'^'>{x) = vr — x). 

Similarly to the case of 90°-walls, the choice of the admissible class An ensures the 
conditions (13) at inhnity, and for other choices of wall orientation there is a net line 
charge as well. 

The remainder of this article is structured as follows. In §3, we present proofs of the 
results given above, using primarily variational methods. In §4, we conduct a numerical 
study of ID domain walls in fourfold materials using ID simulations, and perform 2D 
simulations of a rectangular him element to observe static magnetization conhgurations 
involving these walls. Finally in §5 we conclude and suggest some further extensions to 
this work. 
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3. Proofs of main results 

The following section is devoted to motivating the statements of theorems 1 and 2, and 
presenting their proofs. 


3.1. 90° walls: Proof of theorem 1 

Let us begin by motivating the precise statement of the theorem. Firstly we note that 
in principle, the result one would like to obtain is existence of a solution to (9) which 
satishes (without loss of generality) the conditions 

lim 0{x) = 0, lim 6{x) = 7r/2, (14) 

a:—)-+co x^—oo 

and is in some sense physical, i.e. having hnite energy per unit length of the domain 
wall. Let us recall the energy per unit length in (10) for a ID wall of orientation fi. In 
explicit terms, it reads 


£;.(«) = + 29 ) di 



(sin(6'(a;) — /3) — sm{6{y) — 
(x - vY 


dx dy. 


(15) 


We would like to choose an admissible class of minimizers corresponding to 90° 
transition layers with hnite energy which connect two of the global minima of Sp (given 
by 6{x) = N'k/2 for any G Z). Without loss of generality, we can consider prohles 
satisfying (14). For 9 G iLioj,(M) the local part of the energy is locally well-dehned. 
In order that the nonlocal term be hnite for such prohles, we must constrain the wall 
orientation (3 appropriately so as to avoid incurring a net magnetic charge across the 
wall. To accomplish this we take f3 = —7r/4. 

The ID 90°-wall energy may be expressed as 


= ]^ + dx +^ ||sin(0 + 7r/4)||^i/2(R), (16) 

and the Euler-Lagrange equation associated to £^_ 7 r /4 is now given formally by 

1 / i2 \ 1/2 

0 =+ ^sin46'+ ^ cos(6'+ 7r/4) fj sin(6'+ 7r/4), (17) 


with limit conditions (14). 

The motivation for the statement of theorem 1 should now be clear. We present 
a slightly abbreviated proof of this result; much of the machinery follows directly from 
the work of Chermisi and Muratov [30] concerning Neel walls in uniaxial materials. 
Thus we refer the reader to their work when proving certain steps, and focus here on 
aspects which are signihcantly diherent. 

It is convenient to hrst record some preliminary lemmas. 




One-dimensional domain walls 


Lemma 1 (Restriction of rotations). Let 6 G Then there exists 9 G At^/2 such 

that 0 (R) C [ 0 , 7 r/ 2 ] and S_Tr/4:{9) < i^-7r/4(6*). 

Proof. We let 6 *^ : M —)■ [0, 7 r/ 2 ] be defined, for fc G Z, by 


We then have 



6{x) — kir if 6{x) G [kn, {k + \)'x), 
kir — 9{x) if 9{x) G [{k — |) 7 r, kir). 


sin 2^11 



sin 20 



11^ IIl2(R) - 11 11l2(R) ’ 


and the ineqnality in the lemma follows from setting 6 






= 9^, the ineqnality 


dx'^) 


1/2 


|m| dx, 


and the fact that |sin(0 + 7 r/ 4 )| = sin(0^ + 7 r/ 4 ). □ 

Given this lemma, we may restrict the admissible class to those 9 G At,i 2 also 
satisfying 0(M) C [0, 7 r/ 2 ]. It is then usefnl to define a fnnction p : M —)■ [0, 7 r/ 4 ] G 
corresponding to each 9 in this restricted class, via 


p(x) 


9{x) if 9{x) G [0, 7 r/ 4 ), 

71/2 —9{x) if 0(x) G [vr/d, 7 r/ 2 ]. 


(18) 


It is clear that p satisfies = E-t,iA9). One then has the following lemma: 

Lemma 2 (Coercivity). Let 9 G .4,^/2 he such that 0(R) C [0,7r/2] and let p he defined 
as in (18). Then 

T-^Ia{p) > ^ IIpIIhi(k) + ^ l|sin(p + 7r/4)||pi/2(„). 

Proof. The bound follows immediately from the fact that, since p{x) G [0,7r/4] for all 
a; G R, one has 

cosp(a;) > l/\/2, sinp(x) > p{x)/V2, 

and the identity sin^ 2p = 4 sin^ p cos^ p. □ 

The following lemma provides useful properties for candidate minimizers. 

Lemma 3 (Rearrangement). Let 9 G At^i 2 with 0(R) C [0,7r/2]. Then 39° G 
G(R; [0,7r/2]) satisfying £^_^/4(0°) < £^_^/4(0) and the following properties: 

lim 9°{x) = 0, lim 9°{x) = 0°(O) = j, 9°{x) = ^ - 9°{-x), 

X—)-+oo tc—)•—oo Z ^ Z 


and 9° is non-increasing. 
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Proof. The proof here is similar to that of lemma 4 in [30]. Let p : M —)■ [0,7r/4] 
be dehned as in (18) and let p* denote its symmetric decreasing rearrangement. By 
standard properties of rearrangements (as in [30]) we have hrstly that 


sm 


’ 2p* dx = 


sm 


’ 2p(lx. 


(19) 


Secondly, 

[sin(p + 7r/4) — sin(7r/4)]* = sin(p* + 7r/4) — sin(7r/4). 
Using this, along with lemma 3 in [30], we see that 

||sm(p+ f)||5,i/2(R) > \\Mp* + f)||ii/2(K) • 

From [43, Lemma 7.17], we have 


(PxY < 


pldx. 


( 20 ) 


( 21 ) 


Thus, combining (19), (20), and (21), we obtain 


S-n/A{p*) < S-n/4{p) = S-n/4{0). 
Finally, dehning 6° G C(M; [0,7r/2]) via 


e^x) 


p*{x) if x > 0 

7i/2 — p*{—x) if x < 0, 


( 22 ) 


it is clear that £-n/ 4 { 0 °) = S_T^/ 4 {p*), and that 9° satishes the properties given in the 
lemma. □ 


We now turn to the proof of theorem 1. 

Step 1: Existence. Take a minimizing sequence {6*„} C ^ 7 r/ 2 - By translation- 
invariance and lemmas 1 and 3 we may assume 

G U(M; [0,7r/2]), 6'n(0) = 7r/4, 9n{x) = 7i/2 - 9n{-x), 

and dn are non-increasing. For each n, dehne : M —)■ [0,7r/4] as in (18). From lemma 
2 we have ^ 

^ lIPnll^qR) + ^ ||Mn|lii-i/2(]R) < C < OO, 

where Un = sin(p„-|-7r/4)—sin(7r/4). We may then extract a subsequence (not relabelled) 
such that the following weak convergences hold: 

^ p in 77^ (M), Un ^ u in 77^/^(M), 

Moreover, by compact embedding of the spaces 77^ (M) and 77^/^ (M) into L;^q^(M), we 
have, upon extraction of a further subsequence: 

p„ —)■ p and Un ^ u strongly in LLm and a.e. in R. 
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Therefore, we have u = sin(p + 7 r/ 4 ) — sin( 7 r/ 4 ) a.e. in M. Then, using Fatou’s lemma 
applied to the second term in the dehnition of the energy and lower semicontinuity 
of homogeneous and norms, one obtains weak lower semicontinuity of 

^-wliip) with respect to the weak convergences considered. Passing to the limit n —)■ oo, 
we thus obtain 

^- 7 t/ 4 {p) < liminf£:_^/ 4 (p„) = liminf 

n^oo n—^oo 

Given such a p, we construct a function 6 ^^°^ : R —)■ [ 0 , 7 r/ 2 ] as in (22), that has the 
properties 

0(0)(0) = 7 r/ 4 , 0(°)(a:) = 7i/2 - 0(°)(-x), lim = 0, 

X^+OO 

and 0(0) is non-increasing. Since £)_ 7 ^ 74 ( 0 ) = we then conclude that 0 (o) is a 

minimizer. 

Step 2: Regularity. Since 0 (o) is a minimizer, it must also be a weak solution to 
(17). That is, 'ifj = 0(o) is a weak solution of the equation 


tjj” -\- a{x) cos 2-^ -f- h{x) cos('^ -|- vr/4) = 0, 


(23) 


where a{x) = |sin 20 (o) and h{x) = | sin(0 -|- 7 r/ 4 ). It is easy to see from 

monotonicity of 0(o)(a;) and the limit conditions (14) that a{x) G L^(R). Moreover, 
since 



du 

dx 


2 

1,2 (R) 


Vm e H\R), 


(24) 


and since [sin(0(o) -|- 7 r/4)]' G L^(R), we may see that b{x) G L^(R), as well. Examining 
(23), we may then conclude that 0® G T^(R), and thus that 0® G Furthermore, 

Morrey’s theorem then tells us that 0 i°^ G G(R) fl L°°(R), and so 0(o) g G^(R). Going 
a step further, we differentiate a{x) and b{x) to get a'{x) = sin 20 (°)cos 20 (°) 0 i°^ and 

b'{x) = I f — [(sin(0(°) -I- 7r/4))a,]. Glearly a'{x) G T^(R), so a G 

We would like to show the same for b{x). Applying (24) to u = (sin(0(°) -|- 7 r/ 4 )) 3 ;, 
this then amounts to having to show that G L^(R) (we already know that u G L^(R)). 
We have Ux = (sin(0(°) + 7r/4:))xx = cos(0(°) -|- 7 r/ 4 ) — (0i°^)^ sin(0(°) -|- 7 r/ 4 ). Applying 

the interpolation inequality as in [29]: 



< l|9f 


1/2 

L°°(R) 


^fll 


1/2 

L 2 (R) 


< 00 , 


one may then conclude that b'{x) G L^(R). Again, from (23), this gives 0® G i7^(R) 
and thus 0® G G(R) nL°°(R), implying that 0(°) G G^(R) is a classical solution to (17). 
These arguments may then be bootstrapped to conclude that 0(°) G G°°(R). 

Step 3: Strict monotonicity. The proof here is similar to that in [30]. From the 
previous steps we know that 6^'^ < 0 on R and 0(°) G G°°(R). We claim that in fact 
0® < 0 on R. First, as in [29], we show 6'®(0) < 0: if we assume 0(°)(O) = 7 r /4 
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and = 0, then uniqneness of solntions for the initial valne problem implies 

9^^\x) = 7r/4 identically, contradicting the limit conditions (14). 

Next, we show < 0 on M"*". Assnme there exists an x* > 0 snch that this is false: 
i.e. = 0. This, together with the non-increasing property of 6^^\x) we already 

have, implies that 6^J{x^) = 0 as well. Differentiating the Euler-Lagrange equation and 
evaluating at x*, we then have 


61 ( 0 ) 


= ^cos(6>^°^(x*)-t-7r/4)5((x*), 


.X* 


where we dehned 


9{x) = 


d2 


1 


1/2 


dx^/ 


cos(6*^°^(x) -I- 7r/4). 


Computing (/(x*) directly, using the integral representation of 
that by the assumption on x* the obtained integral converges, we find 


dx2 J 


1/2 


and noticing 


g(x) = - 

TT 


4:X^y9^x\y) cos(6*(°)(i/) -F 7r/4) 
(x* -|/)^(x* +yy 


dy. 


Since 9^^ < 0 on M, with the inequality strict on a set of positive measure, we see that 
g{x^) > 0. Moreover, we know that cos(6*(°) -|- 7r/4) > 0 on M+, and hence 9^xx{x^) > 0, 
implying that 9^^^ is locally increasing at x*. This is a contradiction. Strict monotonicity 
on M follows from the fact that 6'®(x) = 7r/2 — —x) as proved in step 1. 

Step 4- Uniqueness (up to translations). Let 9^1 and 9^'^^ be two distinct minimizers. 
After suitable translations these satisfy = 7r/4. Define u = sin(6'-|-7r/4). 

We may write the energy i^_,r /4 in terms of u as S-Tr/4{9) = E{u), where 

2 . / ^2 \ 1/2 






dx -I— 


u 


\) <^^ + 1 


u 


dx^ / 


u dx 


Then dehne 


9{x) = 


sin ^{\{u^^\x)-\-uU‘\x))) — Ti/A 


37r/4 —sin -|-n^^^(x))) 


,(2)/ 


if X > 0, 
if X < 0, 


where = sin(6'^^i -I- 7r/4) and uU‘'> = sin(6*‘^^i -I- 7r/4). Note that sin(6* -I- 7r/4) = 
and is symmetric decreasing. 

Arguing as in [30], we have 

\^x) S 2 

Then, since the last two terms in E{u) above are strictly convex in u (at least, for 
u G [;^, 1], which is equal to the possible range of values of u that we have here), we 


have that 


S-^/a{9) < 




which contradicts the minimality of 9'dl and 9^‘^\ This concludes the proof of Theorem 

1 . 
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Remark 1. It should be possible to use the arguments of [31] to prove that the 
minimizers in theorem 1 are the unique monotone 90° wall profiles, i.e., that the 
minimizers are the unique, up to translations, monotone critical points of the energy 

7r/4 ^'Kj2' 


3.2. 180° walls: Proof of theorem 2 


Similarly to the above analysis of 90°-walls, we would first like to write down a ID wall 
energy for 180°-walls and define an appropriate admissible class for minimizers. We 
assume, without loss of generality, a prohle connecting the optimal uniform states 6 = n 
at —oo and 0 = 0 at +oo. The appropriate wall orientation to avoid net charge is given 
by /3 = 0. The ID wall energy is thus expressed as 

(|^'r +^sin2 20) dx + ^||sin0|||i/2(R), (25) 


with the admissible class for minimizers given by A^- The Euler-Lagrange equation 
associated to So is given by 


0 


-dxx H— sin46' H— cos 6 

XX 4 2 



1/2 

sin 9, 


(26) 


with limit conditions 

lim 9{x) = 0, lim 9{x) = tt. (27) 

X^+OO X^—OO 

We now turn to the proof of theorem 2. Much of the proof follows by direct analogy 
with the previous section. Indeed, following the proof of theorem 1, it is easy to see that 
lemmas 1 and 3 generalize trivially such that we can immediately restrict the admissible 
class to non-increasing 9 G At^ satisfying 9{R) = (0,7r) along with the properties 

^(0) = 9(x) = TT -9(-x). (28) 

Note that we do not prove uniqueness here; the methods used to prove uniqueness for 
90°-walls in the previous section do not apply due to the fact that the anisotropy energy 
is nonconvex as a function of u = sin 6* for some values of 9 which, in this case, the 
prohle must take. The question of uniqueness of minimizers for this problem remains 
open. 

Lemma 2, however, does not generalize to this case. This is again due to the 
same issue of nonconvexity that causes problems for uniqueness. In order to obtain 
compactness in 77^ (R) for the minimizing sequence, we have to prove a bound on the 
L^(R) norm of p, this time dehned as 


p{x) 


9{x) if 9{x) G [0,7r/2), 

7r — 9{x) \f 9{x) & [tt 


(29) 


such that once again So(9) = So(p). Physically, the issue with compactness which has 
occurred here can be interpreted as the question of whether it is energetically preferable 
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for the 180° transition layer to split into two well-separated 90°-walls (which wonld 
resnlt in the minimizing seqnence weakly converging to 7r/2, which is clearly ontside of 
the admissible class) or whether the fnll transition occnrs mostly over a hnite interval. 

It is clear that the local part of the energy (25) is nnchanged by having an arbitrarily 
large region with 6* = 7r/2. Below, we show that this is not possible dne to the nonlocal 
term (a similar argnment was nsed in the analysis of existence of 360° walls in nniaxial 
materials [40]). Indeed, for 6 in the above class, there exist two nnmbers 0 < a < 6 snch 
that 6{a) = tt/S and 6{b) = n/Q. From the anisotropy term alone, we get that h — a 
remains bonnded above by a mnltiple of the energy. Using the known symmetry of 9 
(i.e. that sin6*(—x) = sin6*(a;)), the nonlocal term in the energy (25) is proportional to 

(sin6*(x) — sin 9{y))'^ (sin6*(x) — sm9{y))'^' 


I = 


[X 


yf 


{x + yy 


dx dy. 


(30) 


We can estimate I from below by neglecting the interval (a, b) from the integrals. 
Dehning 


(sin6*(x) — sin6*(|/))^ (sin 6'(a:) — sin 0(|/))^ 
y) = -7Z— , .A9 - ^ 0) 


(X 


yf 


{x + yy 


(31) 


we have 


n a pa poo poo poo 

f{x,y)dxdy + 2 / / f{x,y)dxdy+ / / f{x,y)dxdy. (32) 

Jo Jb Jb Jb 

We may then estimate the cross term as follows: 


m.y)i.iy>cy jy d. d, = C In . 

for some nniversal C > 0. 

One can see that the nonlocal term forces both a and b to be bonnded by a mnltiple 
of the energy. Then, in order to get a bonnd on ||p|| 2 , 2 (]r), we can nse the following bonnds 
on the anisotropy energy 

POO 2 ^ ^oo 

O > To(p) > T / sin^26*dx = - / sin^ 26*- / sin^26*dx. (33) 

“ 4 io 4 io ^Jb 

For the second term, we have 9{x) G [0,7r/6] for x G [6, cx)), snch that cos6* > \/3/2 and 
sin 6 > 9\/d/2. Using the identity sin^ 26* = 4 cos^ 6* sin^ 9 we then have 



1 

4 



sin^ 29 > 


16 



(34) 


Additionally, 


9^ dx < C. 


(35) 


Jo 

since we know b is hnite and 6*(M) C [0,7r]. We can then conclnde that the transition 
from 0 = 7r/2to6* = Oasa: increases takes place over a hnite distance from x = 0, and 


that 


ll6'llL2(ffi) < C, 


(36) 
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for some C > 0 depending on So{9). 

One then has, for a minimizing sequence (pn), that ||Pn||j:/i(iR) < C and 

||sinp„||^i/ 2 (]g) < O, for some O > 0 and all n G M. One can then extract weakly 
convergent subsequences, and the proof of existence of minimizers, along with the 
regularity and strict monotonicity, follow in precisely the same fashion as in the proof 
of theorem 1. 

4. Numerical study 

In order to assess the role played by the one-dimensional domain wall solutions 
constructed in the preceding sections for the domain patterns in two-dimensional dims, 
we use a hnite-difference scheme to solve the overdamped Landau-Lifshitz-Gilbert 
equation (5), coupled with optimal-grid-based methods to compute the stray held. The 
algorithm is fully discussed in the work of Muratov and Osipov [21], 

4-1. One-dimensional simulations 

We hrst aim to solve (5) in ID to obtain the wall prohles corresponding to our existence 
results in the previous section. We evolve the equation above, beginning from initial 
conditions which approximate the domain wall prohle in question, until a steady state 
is reached. Figure 1 below displays 90° wall prohles for z/ = 1, 5 and 50 (upper panels), 
and their corresponding tails in log-log coordinates (lower panels). Figure 2 displays 
the analogous plots for 180° walls. 

Examining figures 1 and 2, one can observe that the 90° and 180° wall prohles 
behave in much the same manner as z/ is increased. Additionally, the tails all display 
algebraic decay proportional to 1/x^ far from the core, as was proved to be the case 
for 180° Neel walls in uniaxial materials [30]. This decay sets in further away from the 
core as n is increased, and is clearly preceded for large n (e.g. in the plots for z/ = 50) 
by a logarithmic crossover region between the core and the algebraic tails, as is well 
known for Neel walls in uniaxial materials [1,27,28]. We can conclude that these wall 
prohles show ehectively the same behaviour as Neel walls in uniaxial materials. Finally, 
in hgure 2, for u = 1 one can see that the 180° wall is starting to separate into two 90° 
walls. This is due to the fact that at z/ = 0 only the 90° walls exist, while for z^ > 0 the 
180° wall is stabilized purely by the magnetostatic interaction. 

4-2. Two-dimensional simulations 

We now solve (5) on a spatial domain D = [0, L^] x [0, Ly] with the edges of the domain 
aligned with the easy axes of the material. We have ehectively just three parameters 
in the model: L^, Ly and v. Each of these has the same intuitive ehect on the energy: 
increasing domain size or v increases the strength of the magnetostatic interaction, 
relative to anisotropy and exchange. 
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u = I u = 5 u = 50 



XXX 




Figure 1. Computed ID 90° wall profiles for = 1, 5 and 50. Upper panels show the 
wall profiles near the transition layer. Lower panels show the corresponding decay in 
the tails, plotted in log-log coordinates. Red line segments in the lower panels indicate 
an algebraic decay oi 9 ^ 1/x^. 



Figure 2. Computed ID 180° wall profiles for = 1, 5 and 50. Upper panels show the 
wall profiles near the transition layer. Lower panels show the corresponding decay in 
the tails, plotted in log-log coordinates. Red line segments in the lower panels indicate 
an algebraic decay oi 9 ^ \/x^. 
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In the figures below we display stationary configurations (remanent states) for 
a range of domain sizes and values of u. Physically, the fixed parameters we use 
correspond to those of an epitaxial cobalt film, where u can be thought of as the suitably 
rescaled film thickness. The parameters (recall the definitions in equation (3)) are given 
approximately by £ ~ 3.37nm, Q ^ 0.08, and thus L ^ 12nm and u represents the film 
thickness in nanometers [22]. 

Starting from 3 different initial states, we can isolate 4 distinct stationary solutions, 
as follows. In figure 3 below we observe the well known ‘C’ (panel (a)) and ‘S’ (panel (b)) 
states. These states are close to monodomain states, but with edge domains appearing 
along the short edges of the sample to appease the magnetostatic energy at the boundary. 
This energy term prefers the magnetization to align tangent to the boundary wherever 
possible. These states result in magnetically charged regions close to the short edges, 
but with zero energy density in the bulk of the sample. There are necessarily half 
boundary vortices, which also carry charge, in two of the corners. The C state has lower 
energy than the S state. In both panels, the size of the sample is given (in units of L) by 
Lj. = 8,Ly = 16 (the figures are to scale), with z/ = 5. In panel (b), the initial condition 
for the simulation was a monodomain with 6 = 7r/3; in panel (a), we took 6 = 7r/2 in 
the lower half of the sample and 6 = —7r/2 in the upper half. In the figures, the color 
indicates the value of 6{x,y), and the vectors show the corresponding magnetization 
field. 

In figure 4, we see instead states consisting of an arrangement of 90° and 180° 
domain walls, with interesting phenomena appearing at the lower boundary. In both 
figures, one has = 32, Ly = 64, and the simulations were initialized with 6^ = 0 in 
the right half of the sample and 6* = tt in the left half. In panel (a), n = 5; in panel (b), 
z/ = 10. 

In the top halves of both panels, one observes the same phenomenology. There is 
a 180°-wall aligned vertically in the center of sample, which splits into two 90°-walls, 
which then terminate at the top corners. This arrangement separates the top half of the 
sample into three domains with orientations along the easy directions 6 = 0, 7r/2, and tt, 
which also coincide with the orientations of the edges. The arrangements are strongly 
reminiscent of the so-called Landau states [1]; we refer to them as half-Landau states. 

It is simple to understand why the transition from an edge-domain state to a 
half-Landau state is preferred as the relative strength of the magnetostatic interaction 
increases, just by considering the top halves of the figures. The half-Landau states are 
in principle charge free (in the top half), while the edge-domains are not. Thus, while 
sacrificing some exchange energy in order to form the domain walls, the magnetostatic 
energy is reduced. 

In principle as mentioned above, according to the theory of section 3, each of the 
domain walls in the half-Landau state should be individually charge-free. While this is 
true of the 180°-wall in the center, it is not quite true of the 90° walls here. Indeed, 
their orientations are not quite at 45° to the easy axes, and so they are each slightly 
charged (in opposite senses). Since the stray-field energy of such walls would diverge 
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Figure 3. (color online) C-state (panel (a)) and S-state (panel (b)). Domain size in 
both panels is = 8, Ly = 16, with v = 5. 


logarithmically as the size of the domain increased, we would expect that in the large- 
domain limit, they would converge to a 45° orientation. 

Let us now discuss the lower edges of the two half-Landau states pictured in hgure 
4. Firstly, we note that in the model we consider, the magnetization is prevented from 
forming the full Landau state, since in this state there would necessarily be a magnetic 
vortex included in the sample, and these are of inhnite energy in our 2D model. Indeed, 
admissible magnetization states in this model must have topological degree zero (i.e. 
be continuously deformable to the uniform state); a vortex has degree one, as does a 
Landau state. 

In panel (a) of hgure 4, we have u = b. On the lower boundary, the magnetization 
is aligned with 6 = 7r/2—the same as the top boundary. There are boundary vortices 
(quarter vortices) of opposite charges in the lower corners, and what appear to be two 
oppositely charged 90° “boundary domain walls” joining the bulk domains on the left 
and right to the boundary orientation. In panel (b), the magnetization on the lower 
boundary is instead aligned mostly with 6 = — 11 / 2 , antiparallel to the magnetization 
at top boundary, and performs a full 27r rotation in the center of the lower edge (i.e. 
there is a bound pair or boundary vortices there), such that the whole conhguration 
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Figure 4. (color online) Half-Landau states. Domain size in both panels is = 
32, Ly = 64. In panel (a), v = b and the configuration has boundary vortices in the 
lower corners. In panel (b), z/ = 10 and we see a configuration with 2 boundary vortices 
as a bound pair in the center of the lower boundary. 


has degree zero. Additionally, the transition layers at the lower boundary look closer in 
structure to edge domains than to ID boundary walls. 

In panel (a), the charge on the lower boundary is spread out across the whole 
edge while the exchange energy is conhned close to the edge in the boundary walls. 
Conversely in panel (b), the charge on the boundary is focused at the bound pair of 
boundary vortices in the center, while the exchange energy is more spread out across 
the edge domains. For larger i/, it thus appears energetically preferable to concentrate 
the charge in a smaller region. 

5. Conclusions and further work 

We have proven existence results concerning 90° and 180° domain walls, viewed as 
minimizers of the ID domain wall energy, in thin him ferromagnets with fourfold in-plane 
anisotropy. Moreover we are able to learn a lot of information about these structures, 
including strict monotonicity, smoothness, and, in the case of 90°-walls, as well as 
uniqueness. Further problems here include proving uniqueness (or not) of the 180°- 
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wall, and studying walls with more winding such as 360°-walls, or materials with more 
exotic crystalline anisotropy (see, e.g., [44]). We also presented numerically computed 
ID wall prohles corresponding to our existence results, and magnetization conhgurations 
for rectangular samples of fourfold materials, which feature slightly charged walls in the 
bulk and charged walls at the boundary. This investigation poses possible questions 
about domain walls in regimes where the penalty for walls having net charges is relaxed 
such that charged walls, both in the bulk and at the boundary may be allowed in the 
large domain limit, and it would be interesting to study the T-limit of such a 2D thin 
him energy. Additionally, there is the question of existence of ‘boundary walls’ in the 
thin-hlm regime considered in this article: this will be addressed elsewhere [45]. 
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